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ABSTRACT 

This work investigates the short wavelength stability of the magnetopause between a 
rapidly-rotating, supersonic, dense accretion disc and a slowly-rotating low-density magne- 
tosphere of a magnetized star. The magnetopause is a strong shear layer with rapid changes 
in the azimuthal velocity, the density, and the magnetic field over a short radial distance and 
thus the Kelvin-Helmholtz (KH) instability may be important. The plasma dynamics is treated 
using non-relativistic, compressible (isentropic) magnetohydrodynamics. It is necessary to in- 
clude the displacement current in order that plasma wave velocities remain less than the speed 
of light. We focus mainly on the case of a star with an aligned dipole magnetic field so that the 
magnetic field is axial in the disc midplane and perpendicular to the disc flow velocity. How- 
ever, we also give results for cases where the magnetic field is at an arbitrary angle to the flow 
velocity. For the aligned dipole case the magnetopause is most unstable for KH waves propa- 
gating in the azimuthal direction perpendicular to the magnetic field which tends to stabilize 
waves propagating parallel to it. The wave phase velocity is that of the disc matter. A quasi- 
linear theory of the saturation of the instability leads to a wavenumber (k) power spectrum 
oc of the density and temperature fluctuations of the magnetopause, and it gives the mass 
accretion and angular momentum inflow rates across the magnetopause. For self-consistent 
conditions this mass accretion rate will be equal to the disc accretion rate at large distances 
from the magnetopause. 

Key words: accretion, accretion discs — stars: neutron — X-rays: binaries — magnetohy- 
drodynamics — Instabilities — Waves 



1 INTRODUCTION 

This work investigates the short wavelength stability of the inter- 
face or magnetopause between a rapidly rotating accretion disc 
and the slowly-rotating, low-density magnetosphere of a magne- 
tized star. The nature of the magnetopause is sketched in Figure 1. 
The rotating disc matter is "held off" by the star's magnetosphere 
where the magnetic field is strong and the density is small. The disc 
matter rotates at approximately the Keplerian velocity which is typ- 
ically much larger than the velocity of the magnetospheric plasma 
which corotates with the angular velocity of the star. Thus the in- 
terface involves a strong shear layer as sketched in the bottom part 
of Figure 1 . Understanding the instabilities of the magnetopause is 
important for understanding both the transport of matter and angu- 
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lar momentum towards the star and the temporal variability of the 
sources (van der Klis 2006). 

The magnetohydrodynamic (MHD) stability of configurations 
such as in Figure 1 was investigated earlier by Li & Narayan (2003) 
assuming incompressible flow and perturbations independent of z, 
but with no restrictions on the azimuthal wavelength. They found 
both long-wavelength (i.e., < r) Rayleigh-Taylor (RT) and short 
wavelength (<k r) Kelvin-Helmholtz (KH) instabilities for differ- 
ent conditions of the shear layer. Earlier, Burnard, Lea, & Arons 
(1983) studied the short wavelength KH instability of the mag- 
netopause of a spherically accreting rotating magnetized star. Re- 
cently Tsang & Lai (2009) have studied the long wavelength RT 
stability of the sharp interface including the compressibility of the 
media. The MHD stability of the magnetopause for cases where the 
shear layer has appreciable radial width was studied by Lovelace 
& Romanova ((2007) and Lovelace, Turner, & Romanova (2009) 
for compressible, non-barotropic perturbations independent of z but 
no restrictions on the azimuthal wavelength. They found a reso- 
nant long-wavelength Rossby Wave Instability (RWI; Lovelace et 
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Figure 2. Geometry of a small patch of the interface between the star's 
magnetosphere and the external disc flow. 



Figure 1. Sketch of the inner region of a disc around a rotating magnetized 
star suggested by MHD simulations of Romanova et al. (2008). Inside of r„, 
the plasma flows along the magnetic field lines to the surface of the star in 
what is termed a funnel flow (denoted ff in the figure). The bottom part of 
the figure shows the midplane profiles of the the density p, magnetic field 
B z , and azimuthal velocity v*. 



al. 1999) which may contribute to the observed twin kilo-Hertz 
quasi-periodic oscillations of low-mass X-ray binaries (van der Klis 
2006). 

Here we consider a thin magnetopause and short wavelengths 
where the KH instability is expected to be important. There is a vast 
literature on the magnetized KH instability in different astrophysi- 
cal and space applications, to the magnetopause of rotating planets 
in the solar wind (e.g., Miura & Pritchett 1982; Roy-Choudhury 
& Lovelace 1986; Faganello, Califano, & Pegoraro 2008), to the 
stability of astrophysical jets (e.g., Hardee 2007; Osmanov et al. 
2008), to the stability of interfaces of molecular/atomic clouds in 
the interstellar medium (e.g., Hunter, Whitaker, & Lovelace 1998), 
and to the interface of an unmagnetized dense plasma blob falling 
through the strong magnetic field of a neutron star (Arons & Lea 
1980). 

Section 2 of the paper gives the basic equations where the fluid 
motion is assumed non-relativistic but the displacement current is 
retained in Maxwell's equations in order to keep the wave speeds 
less than the light speed. This section also describes the assumed 
equilibrium, the waves in the disc plasma, and the MHD waves in 
the magnetospheric plasma. The disc flow speed is much larger than 
the disc sound speed. Section 3 obtains a dispersion relation for the 
Kelvin-Helmholtz modes and develops an approximates solution 
for cases where the magnetopsheric density (p 2 ) is much less than 
the disc density (pi). Section 4 discusses the non-linear saturation 
of the KH modes and develops a quasi-linear theory for the mass 
accretion rate across the magnetopause. Section 5 discusses briefly 
the case where the magnetic field is at an arbitrary angle relative 
to flow velocity in the disc. Section 5 gives the conclusions of this 
work. 



2 THEORY 

The rotating disc matter is termed region 1, while the magneto- 
sphere matter is termed region 2. In the magnetosphere the mag- 
netic field B is significant whereas in the external medium B is 
negligible as suggested by MHD simulations (Romanova, Kulka- 
rni, & Lovelace 2008; Kulkami & Lovelace 2008). Depending on 
the region, the flow is described by non-relativistic hydrodynamic 
equations or by MHD equations, 

= -Vp+ -JxB+pg, ^+V-(pv) = 0, (1) 



dt 

where g = -V<J? is the gravitational acceleration due to the star 
and O is the gravitational potential. In addition we have Maxwell's 
equations including the displacement current, 

1 <9B 

c~dt 

and Ohm's law for infinite conductivity, 



An 1 <9E 
VxB = — J+ - — 

c c at 



VxE 



(2) 



= E + vxB/c. 



(3) 



where v is the flow velocity, B the magnetic field, p the plasma 
pressure, and V • B = 0. We assume isentropic flow with y = 5/3 
in both media: p = K\p y in the external medium (region 1) and 
P = K iP y in the magnetospheric plasma (region 2). 

One can eliminate E and J from the equations to obtain 

<9B 

~di '' 
and 
d\ 



V x (v x B) . 



(4) 



dt 



V \p + 



8ti 



+ B-V 



J ? Bx| ( vxB) +Pg . (5) 



The less familiar last term of this equation arises from the displace- 
ment current in equation (2). It is required in order to have the wave 
speeds less than the velocity of light (Jackson 1975). 

2.1 Equilibrium Flow 

We consider an axisymmetric time-independent equilibrium 
plasma flow. The flow velocity v = v^(r)0 = rQ.,p(r)<f> . That is, 
the accretion velocity v r and the vertical velocity ;( ; are assumed 
negligible compared with v^. Initially, we use an inertial cylindri- 
cal (r,(f>,z) coordinate system. The equilibrium magnetic field is 
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Figure 3. Phase velocities u+ = u ± /|K| (fast and slow) of sound waves 
in the external medium viewed from a reference frame where the magne- 
tosphere is at rest. For the case shown the external medium moves with 
velocity v x = 0.01 and the sound speed is c\ = 0.003 with speeds in units 
of c. The fast/slow phase velocities go through the ;t-axis at v x ± ci. 

B = B(r)z. In §5 we discuss the general case where the magnetic 
field has a general orientation with respect to the flow velocity in 
the disc. The equilibrium flow satisfies -pg = —d[p + B 2 /(8n)]/dr. 
Here, p is the density, p the pressure, and Q. K the Keplerian angular 
rotation rate of a single particle. We assume a pseudo-Newtonian 
potential where S1 2 K = GM/[r(r - r s ) 2 ] with r s = 2GM/c 2 the 
Schwarzschild radius of the star. The effective radial gravitational 
acceleration is g = -r(£l 2 K - H?). 

We consider the case where the plasma properties undergo a 
rapid change at the magnetopause radius r m as sketched in the bot- 
tom of Figure 1. The change in values is assumed to occur over a 
distance Ar„, <k r m . Inside of r m the magnetic field is B, the den- 
sity is p2, the sound speed is ci, and the flow velocity is rQ.,(p , 
where Q., is the angular rotation rate of the star. Outside of r„, the 
magnetic field is negligible, the density is pi, the sound speed is ci, 
and the flow velocity is close to Keplerian, v K (f> . We consider the 
case suggested by 3D MHD simulations (Romanova, Kulkarni, & 
Lovelace 2008; Kulkarni & Romanova 2008) where pi <k pi and 
C2 ~ c\. For such conditions the force balance at r m requires that 
B 2 /$n* Pl =p l c]/j. 

We focus our attention on the stability of short wavelengths 
X = A/2n <k r. This allows us to consider the stability of a 
small patch of the magnetopause separating the two plasmas as 
indicated in Figure 2. The change of coordinates, from cylindri- 
cal to Cartesian, is (r — » z, rip — > x, z — > y). For the considered 
short wavelengths the gravitational force can be neglected in equa- 
tion (5). That is, we consider the Kelvin-Helmholtz rather than the 
Rayleigh-Taylor modes. We choose a reference frame moving with 
the magnetospheric plasma. Consequently the external disc plasma 
moves with a velocity v = v v x with v x = (v K - r£l,) r=rm . 

2.2 Waves in the External Medium 

In region 1 there is no magnetic field so that the small amplitude 
waves ~ exp(/K • r - iojt) in the rest frame of the medium con- 
sist of both sound waves and non-propagating disturbances such as 
variations in the vorticity (Landau & Lifshitz 1987). The dispersion 
relation for the sound waves is oj 2 = K 2 c 2 , where u) is the angular 




Figure 4. Phase velocities u of slow magnetosonic (sms) and fast magne- 
tosonic waves (fms) in the magnetosphere. The magnetic field direction is 
assumed to be along the v-axis. In the top panel of the figure has Va = 1 
and c„, = 0.2, while in the bottom has va = 0.2 and c m = 0.2, where the 
speeds are in units of c. 

frequency, K = (K x , K y , K z ) is the wavevector, c\ = (dp/dp)^ 2 is 
the sound speed in the external medium. In the reference frame we 
use, the external medium is moving with uniform velocity v, so that 
the dispersion relation in this reference frame is 

(oj - K x v x f = K 2 c] . (6) 

which gives oj ± = K x v x ± |K|c sl . For the non-propagating distur- 
bances, cj = K x v x . Figure 3 shows a polar plot of the phase veloc- 
ities u+ = tj ± /\K\ > and uq = a>o/\K\ of these waves as a function 
of (f> which is the angle between K and the x-axis. 

2.3 MHD Modes in the Magnetosphere 

The dispersion relation for the two compressible MHD modes in- 
cluding the displacement current is 

oj 4 - co 2 [K 2 (c 2 + v 2 A ) + (K ■ y A ?c 2 Jc 2 ] + K 2 c 2 (K ■ v A ) 2 = , (7) 
where 

5 » - 1 — ~^7T ' where Cs = ' (8) 

l+v-/c 2 \dp) s 

(e.g., Jackson 1975). Here, oj is the angular frequency of the wave 
and K is the wavevector. Note that the bold-faced \ A 's in equation 
(9) do not have over-tildes and can be larger than c. 

The magnetic field in the magnetosphere is strong and the 
plasma density low so that we expect to have » c 2 or equiva- 
lently v\ » c 2 . In this limit the slow magnetosonic wave branch of 
equation (7) has u) 2 m ss K 2 c 2 cos 2 (^>), where ip is the angle between 
K and B. For this wave the velocity perturbations are parallel to B 
with the result that the wave speed is independent of B. The other, 
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Figure 5. Phase velocities u of the Alfven wave in the magnetosphere as 
well as the fast and slow sound waves u+ in the external plasma. In this 
figure, Cji = c S 2 = 0.0391, \'a = 0.805, and v v = 0.167 all in units of c. The 
dashed circle is v x sin(y) (for < ip < n) and the dashed lines delimit the 
wavevector region of the Kelvin-Helmholtz instability discussed in §3.2. 



fast magnetosonic branch has w 2 m « K v A with u., r 



1 U} m » <4m- This 

wave has velocity perturbations in the plane formed by K and B. 
Both waves give density variations. There is also a non-propagating 
wave analogous to that in region 1 . Figure 4 shows a polar plot of 
the phase velocities of the waves obtained from equation (7). 

The dispersion relation for the shear Alfven wave including 
the displacement current is 



(K-v A ) 2 



where 



(1+V 2 / C 2 )'/ 2 ' 



(9) 



(10) 



where \a = B/(47rp)'' 2 is the usual, non-relativistic Alfven veloc- 
ity which may be larger than c (e.g., Jackson 1975). The velocity 
perturbation of this wave and the field perturbation <5B are perpen- 
dicular to both K and B. Thus the wave does not change the fluid 
density. 



3 LINEARIZED EQUATIONS FOR THE KH MODES 

The geometry of a small patch of the interface between the external 
medium and the magnetosphere is shown in Figure 2. The interface 
is in pressure equilbrium so that p + B 2 /8w is continuous across it. 

Perturbations of all scalar quantities have the form 
/(z) exp(ik • r - icot) for all z, where r = (x, y, 0), k = (k x , k y , 0) is 
the wavevector in the plane of the interface, and the z-dependence 
remains to be determined. The linearized continuity equation is 



iAcoSp + ipk ■ 6v + — (p6v z ) 
dz 



0. 



(11) 



wave frequency in the frame comoving with the external medium. 
The step-function dependences of p(z), v x (z), B x (z), and B y (z), gives 
rise to various delta-function terms in the linearized equations. For 
example, from the top side (region 1 , z > 0) the vertical displace- 
ment of the interface is fei and the velocity at z = e + is 



Sv, 



d d\ 



-iAcoSzi 



while from the bottom side (region 2) it is 

Sv-o = ^-Szi = -iu)5z2 ■ 
at 

We must have dzi = Szi so that 

Sv z {z = e + ) _ Sv z (z = e~) ^ 
to - k x v x (e + ) to 

(Chandrasekhar 1961), where e ± denotes an arbitrarily small posi- 
tive or negative quantity. Thus Sv- is discontinuous across the inter- 
face. Considering perturbations proportional to exp(ik : i^z), equa- 
tion (11) can be written as 

p(z)K-<5v iSv z dp 



dp 



to ~ k x v x (z) cl> - k x v x (z) dz 



(13) 



Here, K = (k x , k y , k z ) is the three-dimensional wavenumber which 
also comes up later. The second term of this equation gives a delta 
function dependence, —6zi(pi — P2)8(z)- 

Linearization of the Euler equation (5) gives 

SB, B • <5v 

- iWpSv x = -ik x Sp + i(k ■ B)- icoB x — — — 

An 47TC 1 



dv x 1 dB x 

+iA(±>p 6z-i -SB- , 

dz An dz 

SB,, B ■ Sv 

iWpSv, = -ik v Sp + f(k • B) — ■ — itoB v — 

An 4nc 2 

1 dB x 

+ ——^SB Z , 
An dz 

d SB- 
iWpSv- = - — Sp + i(k-B)- 



dz 



An 



W = Alj + a>- 



(14) 

(15) 
(16) 



where 



and 



Sp = Sp + —(BySBy + BySBy) . 

An 

The terms oc 1/c 2 are due to the displacement current and they 
must in general be retained. Recall that B = B x x + B y y so that 
K B = k B 

Linearization of the induction equation (4) for the magnetic 
field gives 

Sp k B k • B dv x dB x B K dp 

SB X = B x — - —6v x - -i& - — -Sz + — -j-Sz , 

p Acj Aw dz dz p dz 

SB V = S, — dVy — Sz + -Sz, 

p Act) dz p dz 



SB Z = f'(k • B)Sz ■■ 



k B 

All) 



-5v, 



(17) 
(18) 
(19) 



where <5v = (Sv x , Sv y , 0) and Aco = to - k,v x is the Doppler shifted 



At the interface the terms involving d(v x ,p,B x ,B y )/dz give delta 
functions ~ <5(z). Away from the interface these terms vanish. 

Away from the interface we can combine equations (14) and 
(15) and (17) and (18) to obtain 
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pB .<5v= ^(k.B)(<5p-vl<5p), 
where 



(20) 



F = Aco 2 - (k • y A f , 

and \ A = B/ -y/4/rp. Note K • S\ = Aco(Sp/p) and that K • SB = 0. 

As mentioned we use the reference frame comoving with the 
magnetosphere so that in region 2 Aco = co. From equations (14) 
and (15) we have 



co + \ pK -Sv = K 2 Sp + — 2 (k B)B <5v. 



From equations (17) and (18) we have 
B -S\. 



7 Sp k B 
B« = B — 

p Aid 



(21) 



(22) 



Using the above relation between Sp and dp, and the relation for 
isentropic perturbations Sp = c 2 Sp, we obtain 



Sp- 



(k • y A ) 2 c 2 s 
Aco 2 



Sp. 



(23) 



Equations (21) - (23) can be readily combined to give the dispersion 
relation for the magnetospheric modes. 

Away from the interface equations (11) and (15) give 



<P_ 
dz 2 



Sp = -k 2 z Sp . 



(24) 



For the external disc medium, z > 0, we must have Sp oc exp(ik zl z) 
with 3(k zl ) > so that the perturbation decays as z increases. In 
this region we have 



-k 2 + 



(co - k x v x ) 2 

c 2 
Si 



(25) 



which is identical to equation (6) with K 2 = k 2 + k\ as it should be. 
Dividing this equation by k 2 and taking the square root gives 



k z i(u) 



' [u - sin( V K] 2 



1/2 



(26) 



where k z = k z /\k\, ip is the angle between k and B, and u = w/|k| 
is the phase velocity of the perturbation which is complex with a 
positive imaginary part for an unstable perturbation. The choice of 
the sign in this expression is determined by the condition 3(fci > 0. 

For the magnetospheric plasma, z < 0, we must have Sp oc 
exp(ikaz) with 3(k z2 ) < so that the perturbation decays as -z 
increases. We find 



e = _ k2 + "> 2 -(k-v A )V 2 / c 2 ] 



co 2 c) - (k • y A ) 2 c 2 , 



(27) 



where v A is given by equation (8), c s is given by equation (10), 
and c 2 = c 2 2 + v 2 A is the fast magnetosonic wave speed. Note that 
equation (27) with K 2 = k 2 + k\ is identical to equation (9) as it 
should be. Dividing this equation by k 2 and taking the square root 



gives 



k z2 (u) = : 



u 2 \u 2 - [cos(^)v A ] 2 c 2 2 /c 2 



,1/2 



u 2 c f - [cos.(tp)v A ] 2 c 



- 1 



(28) 



The choice of the sign in this expression is determined by the con- 
dition 3(k z2 ) < 0. 



3.1 Fundamental Dispersion Relation 

Equation (16) can be rewritten as 

^6p = -[wAco-(k-v A ) 2 ]p5z, 



(29) 



where the terms on the right-hand-side are finite. Thus we have 
Sp(z = e + ) = 5p(z = e~). Using the fact that Spu oc exp((fc z i j2 z), 
we find ik z2 Sp = -[WAco-(k ■ \ A ) 2 ]p 2 Sz 2 and ik zl Sp = -Aco 2 pi$zi- 
Taking the ratio of these equations gives k z2 /k zi = (p 2 /pi)[co 2 (l + 
v 2 A /c 2 ) - (k • y A ) 2 ]/(Aco) 2 . This gives 

(30) 



Pi 



V-fc,v,) 2 = g(l + j)p- ( k.v A ) 2 ]. 



This is the fundamental dispersion relation. It agrees with the result 
of HWL in the limit c — » oo. Equation (30) can be rewritten in terms 
of u — co/\k\ as 

F(u) = [u-sm(cp)v x ] 2 -g 2 j^- [u 2 - [cos(^)v A ] 2 ] = , (31) 
[k z2 (u) J 

where k zl (u) is given by equation (26) and k z2 (u) by equation (28), 
and 



B 2 



P2 



pi 4^ Pl c 2 



(32) 



The choice of signs for k z \ (equation 26) and k z2 (equation 28) is 
fixed by the above mentioned requirements that 3(£ z i) > and 
3(£ 2 ) < 0. 

The pressure balance across the interface gives 



Bj; p 2 c\ 

87T y 



(33) 



where c 2 is the sound speed in the magnetosphere. This is the same 
as 



2 2t 2 



2c 



_ _ fh^i 
7 P2 J 
Hence equation (32) can be written as 

2c? = P2 
Pi ' 



pi \ yc ) yc 



(34) 



(35) 



where the last equality takes into account that the plasma motion is 
assumed to be non-relativistic. 



3.2 Approximate Instability Solution for g 2 <sc 1 

We are interested in conditions where j 2 « 1 where an approx- 
imate solution to equation (31) can be developed as follows. The 
two terms of equation (31) are written as 

F = F (u) + g 2 F 1 (u) = 0. (36) 

We develop a perturbation expansion for this equation based on the 
small parameter g 2 . Thus, we take u = u + Su with u is chosen to 
give Fq(uq) = and \Su\ <k \uq\ assumed. Hence 

« F (uq + Su) + g 2 Fi(u + Su) , 
dF 1 d 2 F 

« F («o) + -j±6u + ^—^(Su) 2 + g 2 Fdu ) + 0(g 2 \Su\) . (37) 
du 2 du 2 

We choose uq such that Fq(uq) = which implies that uq = 
s.in(tp)v x . Then we have dF /du = and d 2 F /du 2 ) = 2. Equation 
(37) then gives 
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Su = ±g 



hj(u ) 



k z2 (u ) 



1/2 



[ sin 2 (y)v 2 - cos 2 ((p) 



vl\ 



(38) 



For u = sin((p)v x , equation (26) gives k z i(u Q ) = ±i. Also, for 
u = s'm(<p)v x and for v 2 » c 2 2 and v 2 <K cj, equation (28) gives 
ka( u o) = +i- The mentioned conditions on the imaginary parts of 
the k z 's then implies that k z \lk Z 2 = - 1 • Therefore, equation (38) im- 
plies instability for v x > v A {lwnpY l which is the condition for the 
Kelvin-Helmholtz instability for the considered equilibrium. The 
real part %(u) = %(u>)/k corresponds to the x-component of the 
phase velocity of the perturbation matching flow speed of the ex- 
ternal medium v x . 

For w = sin(</?)v r and Su given by equation (38), one has in 
region 1, pdv zl = k x Spi/a>i from equation (14). From the equa- 
tion following equation (11), we have 6z\ = <5v z i(e + )/w,-, where 6z\ 
is the displacement of the interface. Consequently, 5p\(e + )lp\ = 
oj 2 Sz\ I {k x c 2 sl ), which shows that the perturbation in region 1 in- 
volves a change in the density. The perturbation is a sound wave 
evanescent in the z-direction with amplitude cc exp(ik ■ x - \k z \z - 
ia) r t) where a> r = k x v x and \k z \ = \k x \[l + (gv x /c si ) 2 ] [/2 from equa- 
tion (25). The instability results from the interaction of this wave 
with the Alfven wave (§2.3) in region 2 which has Sp2 = and 
does not change the magnitude of the magnetic field. Perturbations 
which change the magnitude of B are suppressed because they in- 
crease the magnetic energy of the system. For this reason the slow 
and fast magnetosonic waves are not excited. 

Having the wavevector along the magnetic field corresponds 
to bending the field line which requires energy and is stabilizing. 
Thus the maximum growth rate occurs for if = 90° (where k = \k x \) 
with no field line bending and is 



max(wj) = 3(kSu) = gkv x 



1/2 



kv x 



(39) 



This formula for cjj applies only for a restricted range of k: for 
say k > 3r m l where the planar description of the interface is valid, 
and for k < (Ar m ) _1 , where Ar m is the radial thickness of the in- 
terface. The maximum growth rate does not depend explicitly on 
the value of the magnetic field. However, it depends implicitly on 
B since the field allows conditions with p 2 /pi ■« 1 and the field 
enters the expression for k^iu). The same formula for the growth 
rate is found by Li & Narayan (2004) who assume incompressible 
fluid motion in both media. In our treatment the response in the 
low-density magnetized region 2 is an Alfven wave which is in- 
compressible, but the response in the high-density unmagnetized 
region 1 is incompressible only in the zeroth approximation where 
g = in equation (36) and Fq(u) = 0. In the first approximation in- 
cluding the g 2 term the medium is compressible. As discussed in §4 
the compressibility of the region 1 medium gives rise to observable 
fluctuations in the emissions from the interface. 

For comparison with the results of Arons & Lea (1980) who 
took into account the displacement current (important for v A /c > 
1), we need to give the correspondence of our variables with theirs. 
We find that our P2/P1 corresponds with their p e /pt <K 1, that our 
v A is the same as theirs, v A = \B\(4np e y l/2 = \B\(4np 2 y 1 ' 2 % 1, and 
that their a = v A (p e /p b ) [l2 <K v A . In their equation (A48) we can 
neglect terms involving p e /pb in comparison with unity. In the limit 
where a 2 « 1, their equation gives our equation (38) multiplied by 
(1 + v\/c 2 ) 112 . Thus the two results agree only in the limit v 2 /c 2 <<c 
1. For a > 1 the two expressions are also different. 

As a numerical example consider c s i = C& = 0.0391, v A = 
1.36, v A = 0.805, v x = 0.167 (all in units of c), g 2 = p 2 /pi = 
0.00283 and g = 0.0532, Mach number Mi = v x /c si = 4.26, where 



we find (Of = 0.00887fcc for propagation in the x-direction. There is 
instability for <p > arctan(v A /v J ) = 78.3° up to <p = 90°. The range 
propagation directions of the unstable waves is shown in Figure 5 as 
marked KH. The saturation of the exponential growth is discussed 
in §4. 

We can express equation (39) as w; = (pil p\) x,2 k x c s \M\ 
where M\ is the Mach number of the external plasma flow. In 
the commonly considered case of Kelvin-Helmholtz instability be- 
tween equal density media there is instability for waves propagat- 
ing parallel to the flow only for M < 2 3/2 (e.g., Hunter & Whitaker 
1989). In the present case there is no similar limit on the Mach 
number for the growth of waves propagating parallel to the flow. 
We of course have the limit v 2 <c c 2 because of our assumption of 
non-relativistic fluid motion. 



4 NONLINEAR EFFECT OF UNSTABLE KH MODES 

A large number of computer simulation studies have been done on 
the nonlinear evolution of the KH instability for different initial 
configurations (e.g., Zhang, MacFadyen, & Wang 2009; Keppens et 
al. 1999; Frank et al. 1996). The studies do not address the configu- 
ration considered here, but they suggest the approximate treatment 
discussed below. A comparison is made with the nonlinear model 
of Burnard et al. (1983). 

Starting from a sharp interface, exponential growth of the KH 
ceases at some time when typically one or more "Kelvin's Cat 
Eyes" form. Subsequently the interiors of the Eyes may become 
highly irregular but their widths remain roughly constant. For sim- 
plicity, we consider the region of the magnetopause outside of r m 
or equivalently z > (region 1). Also we assume that the wave 
propagates in the x-direction where the growth rate is a maximum. 
From §3 it is clear that the unstable wave initially corrugates the 
interface displacing it by the amount 6z = iSv z [(e + )/ Aa>. The expo- 
nential growth will cease when the width of the Cat Eye is of the 
order of the reduced wavelength. That is, |<5z| = qX, where q =const 
is a number of the order unity and A = l/k. At later times the ini- 
tially sharp interface is effectively smoothed over a distance ~ Sz 
giving the interface a radial thickness Ar m = Sz (Frank et al. 1996). 
This finite thickness acts to stabilize the KH instability for a given 
k = A' 1 . At saturation |<5v z i| = qujk since Aa) = From equa- 
tion (14) we have p\Sv z = i\k\Sp/Aa>. This gives Sp/pi = q(ti)j/k) 2 
and Sp/pi = qw 2 l(kc s i) 2 = qg 2 (v x /c si ) 2 at saturation. The use of 
the linear relations between the fluid variables is plausible in the 
nonlinear regime if \Sp\lpi = qg 2 (y x /c sl ) 2 1. 

In general, there is a spectrum of saturated KH waves with the 
contribution to \Sp\ 2 from wavenumbers of order k equal to k\Spi,\ 2 , 
where |<5p^| 2 is the wavenumber power spectrum of Sp. From the 



previous paragraph we have 



\Spi 



2 2 4 

,2 p\q g 



t'.si 



(40) 



The total mean-squared density fluctuation is therefore (|(5p| 2 > 



J dk\Sp k 



P 2 q g (Vx/c.si) hi(fcmaxAmin)- Here, k mm is In over 



the longest wavelength which is of the order of the size of the 
patch of the interface (<k r m ), and £ max is 2n over the short- 
est wavelength which may be the radial thickness of the magne- 
topause discussed further below. Equation (40) corresponds to a 
"pink-noise" or "one-over-/" spectrum. For the assumed isentropic 
equation of state there are corresponding temperature fluctuations 
\$T k \ 2 /T 2 = (y - l) 2 \Sp k \ 2 lp 2 which can give spatial variations in 
the radiation from the optically thin regions of the magnetopause. 
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Figure 6. Sample realization of the saturated interface for q = 1 and 
^max/^min = 50 viewed in the reference frame comoving with the disc 
plasma (p\). The Fourier amplitudes Spk are generated with a random phase 
and magnitude with rms magnitude oc k^ 2 . 



Connecting the spectrum of waves with observed noise spectra of 
X-ray sources (van der Klis 2006) is, however, beyond the scope of 
this work since it involves integration over the entire magnetopause 
as well a treatment of the radiation transfer. 

Figure 6 shows a sample realization of the interface with the 
wavenumber spectrum given by Equation (40) in a reference frame 
comoving with the disc plasma. 

The equilibrium flow of §2.1 does not include accretion. How- 
ever, the unstable KH modes of §3.2 can give rise to accretion in 
the vicinity of the magnetopause radius r m . The mass flux per unit 
(x, y) area (A) in the z-direction across the magnetopause is simply 



SMkh 
SA 



= %{SpSv" z ) , 



(41) 



where dp = p(k x Sv z + k z Sv z )/Aaj from equation (13), pSv x = 
k x Sp/Aw, and pSv. = i\k z \6p/ Aa> from equations (14). Substitution 
gives 



SM K 



2k 3 c 2 j\Sp\ > 



SA 



2q gp\v x \n 



(42) 



where we have used equations (23) and (39). 

Equation (42) corresponds to an effective accretion speed 
at the magnetopause of KjjSjJ = 2q 2 gv x ln(£ max /£; min ). The accre- 
tion speed proposed by Burnard et al. (1983) agrees with if 
q = 0.11 for k max /k min = 50. For the values given below equation 
(39) and q = 0.1, the accretion speed is 0.00069c. For comparison, 
in a standard alpha-disc model (Shakura & Sunyaev 1973) the ac- 
cretion speed is of the order of w disc = a(c s /vg) 2 Vx, where v K os 
the Keplerian velocity. We find as 10 -3 c assuming a = 0.1, 
c s = Cji, and v K = v x . On the other hand for conditions where 
q = 0(1) and «j£ > «disc> we suggest that the profiles of the den- 
sity and other variables (shown in Figure 1) become less steep hav- 
ing a width Ar m . The present treatment of the KH instability re- 
mains valid for reduced wavelengths X > Ar,„ which corresponds 
to k < k max = (Ar,„)~'. Clearly, for less steep profiles the factor 
ln(£ ra ax)/£min) decreases as does u^. Thus self-consistent condi- 
tions can evolve naturally to give = M disc . 

The unstable KH waves also gives a radial influx of angular 



momentum (about the z-axis) per unit area of the magnetospheric 
patch. For simplicity we neglect the rotation of the magnetosphere, 
(Q,,r m ) 2 <K v 2 K (r m ). Then, 

SF KH 



SA 



r m [%{6pv x Sv* z ) + %(pSv x 6vl)] . 
SMkh 



SA 



(43) 



where the term %(pSv x Sv z ) vanishes and where r m v x is the specific 
angular momentum of the external disc matter outside the magne- 
topause. 



5 GENERAL ORIENTATION OF THE MAGNETIC 
FIELD 

For the general case the star's magnetic moment fi is not aligned 
with fl,. However, we assume the rotation axis of the disc is 
aligned with 17,. The shear layer between the disc and star will 
be inherently time dependent due to the star's rotation. For exam- 
ple, for an orthogonal rotator where fi is in the equatorial plane, the 
plasma in the disc sees a magnetic field reversing its direction with 
an angular frequency a> B = v x /r„, = (v K - rQ.,) r=rm . This frequency 
may be larger or smaller than the growth rate of the KH instability. 
From equation (39) we have a>i/o)B = g(r m /X). For CJi/a}B > 1 sig- 
nificant wave growth can occur before the moves from a given field 
region. In the other limit the wave growth may be recurrent as a 
given plasma region returns to a given field region. In the following 
we consider B to be time-independent. 

For a general field orientation, B is still in the plane of the 
shear layer, that is, the (x, y) plane of Figure 2. An equilibrium with 
a B z component is not possible. Note however that the z-direction 
is not necessarily in the f-direction. For this case we discuss the 
needed modification of Figure 5 where it is appropriate to use the 
angle <p rather than <p. The flow velocity remains in the x-direction 
so that the fast and slow sound wave curves remain unchanged as 
does the dashed circle. The dashed circle can be written as v x cos(</>) 
for -7r/2 < tp < <p. What changes in Figure 5 is that the "fig- 
ure eight" curve for the Alfven wave is rotated by angle 6 say in 
the counter-clockwise direction. The figure-eight curve is give by 
V4I sin(0 - S)\. In place of equation (38) we find 
n 1/2 



Su = ±g 



k zl (u ) 



cos {<t>)v\ 



■ sin 2 (0 



(44) 



k z2 (u )_ 

The limits on the directions of the unstable Kelvin-Helmholtz 
modes are given again by the intersection of the dashed circle in 
Figure 5 and the rotated figure-eight curve. These limits are easily 
found to be 



61 2 = arctan [ tan(0) ± 

\ v A |cos(0)| 



(45) 



Figure 7 shows the dependence of the real and imaginary parts of 
the KH most unstable mode on the tilt angle of the field 9. 



6 CONCLUSIONS 

This work investigated the short wavelength (A <K 2nr m ) stability 
of the magnetopause at radius r m between a rapidly-rotating, super- 
sonic, dense (p 2 ) accretion disc and a slowly-rotating low-density 
magnetosphere (pi » pi) of a magnetized star. The magnetopause 
is a strong shear layer with rapid changes in the azimuthal veloc- 
ity, the density, and the magnetic field over a short radial distance 
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Figure 7. Dependences of the real part of the phase velocity uq = 0) r /k 
and the imaginary part 3{<5«) = cjj/k (in units of c) on the tilt angle of the 
magnetic 6 for the most unstable wave. The parameters are the same as for 
Figure 5. 

(Ar m <K r m ), and thus the Kelvin-Helmholtz (KH) instability may 
be important. This work has focused on the case of a star with an 
aligned dipole magnetic field so that the magnetic field is axial in 
the disc midplane and perpendicular to the disc flow velocity. For 
the aligned dipole case the magnetopause is most unstable for KH 
waves propagating in the azimuthal direction perpendicular to the 
magnetic field. Propagation not perpendicular to the magnetic field 
changes the magnetic field which gives a stabilizing effect. The 
growth rate of the instability is g>; = {pil pi) lll k x v x , where v x is the 
velocity of the disc plasma relative to that of the magnetospheric 
plasma which corotates with the star. The wave phase velocity is 
that of the disc matter. 

We discussed the non-linear saturation of the instability which 
we argued occurs for a mode of wavelength A = 2n/k when the dis- 
placement of the interface between the disc and magnetospheric 
plasmas is of the order of A/2n. From this we developed a quasi- 
linear model which led to a wavenumber power spectrum oc Ar 1 
of the density and temperature fluctuations of the magnetopause. 
The quasi-linear model gave the mass accretion and angular mo- 
mentum inflow rates across the magnetopause. The mass accre- 
tion rate (per unit area of the magnetopause) was found to be 
2(P2/Pi) I/2 Pi v '.vln(fc m axAmin), where £ m a X , mi „ are the maximum and 
minimum wavenumbers discussed in §4. For self-consistent condi- 
tions this mass accretion rate will be equal to the disc accretion rate 
at large distances from the magnetopause. 

We also considered the case where the magnetic field is not 
perpendicular to the flow velocity but tilted by an angle 9 relative 
to the z-axis, where = corresponds to the aligned rotator case 
treated earlier. We found that the maximum growth rate and the 
associated wave phase velocity decrease monotonically with 6 and 
that both approach zero as 6 — > 90°. That the growth rate goes 
to zero is a result of the stabilizing effect of the magnetic field. 
Thus an orthogonal rotator where the star's magnetic moment fi 
is perpendicular to the star's rotation axis fl, is stable to the KH 
mode. 
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